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ABSTRACT 

Accurate stellar population synthesis models are vital in understanding the properties and formation 
histories of galaxies. In order to calibrate and test the reliability of these models, they are often 
compared with observations of star clusters. However, relatively little work has compared these 
models in the ugriz filters, despite the recent widespread use of this filter set. In this paper, we 
compare the integrated colors of globular clusters in the Sloan Digital Sky Survey (SDSS) with those 
predicted from commonly used simple stellar population (SSP) models. The colors are based on 
SDSS observations of M31's clusters and provide the largest population of star clusters with accurate 
photometry available from the survey. As such, it is a unique sample with which to compare SSP 
models with SDSS observations. From this work, we identify a significant offset between the SSP 
models and the clusters' g-r colors, with the models predicting colors which are too red by g-r'^0.1. 
This finding is consistent with previous observations of luminous red galaxies in the SDSS, which 
show a similar discrepancy. The identification of this offset in globular clusters suggests that it is 
very unlikely to be due to a minority population of young stars. The recently updated SSP model 
of Maraston & Stromback better represents the observed g-r colors. This model is based on the 
empirical MILES stellar library, rather than theoretical libraries, suggesting an explanation for the 
g-r discrepancy. 

Subject headings: globular clusters: general - galaxies: star clusters: individual (M31) - galaxies: 
stellar content 



1. INTRODUCTION 

Stellar population synthesis (SPS) models are an im- 
portant astrophysical tool for interpreting photometric 
and spectroscopic observations of galaxies. Modern sur- 
veys can detect vast numbers of galaxies across a range of 
redshifts and can provide important observational con- 
straints on star formation histories and galaxy evolution. 
However, it is generally not possible to resolve the stel- 
lar populations of these galaxies. Instead, their proper- 
ties are often estimated by comparing their integrated 
emission with that predicted from composite stellar pop- 
ulation models. Such models combine our knowledge of 
stellar evolution with stellar spectral libraries to predict 
the integrated emission of an unresolved stellar popula- 
tion. In their simplest form, simple stellar population 
(SSP) models, SPS models predict the integrated emis- 
sion from a stellar population with a single metallicity 
and star formation epoch. The ability of SPS models to 
accurately represent real stellar populations has a direct 
influence on galactic studies. 

For many years, star clusters have provided some of 
the best locations to calibrate parameters in SPS mod- 
els (not all of which can be derived theoretically) and 
to test their ability to reproduce the observed proper- 
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ties. This is because the relatively simple stellar popula- 
tions of these clusters can be directly compared with such 
models. Photometric and spectroscopic observations of 
the Milky Way's (old) globular clusters, (young) open 
clusters and the Large/Small Magellanic Cloud's clusters 
have previously been used to calibrate and test suc h mod- 
els (e .g. Renzini & Fusi Pecci 1988; VandcnBerg fc Clem 
[2OO3I; Mar astoni ,2005^: .An et al.. .2008; .Lvubenova et al. 
2010; Riff el et al.l I20m These studies, which often 
consider observations obtained in the Johnson- Cousins 
UBVRcIc filters, have shown reasonable agreement be- 
tween models and observations. Howev er, se veral dis- 
crepancies, most notably in the B-V (e.g. IWorthev„.199'l: 
MarastonI I2005f) and near- infrared (e.g. ICohen et ail 



2007t IConrov fc GunnI l20Toh colors, have been demon- 



strated. 

Despite the increasing popularity and use of the SDSS's 
ugriz filter set (Fuku gita et al.l 119961) , relatively little 
work exists which compares SSP models to star clusters 
in these filters. This is due primarily to a relative lack 
of cluster photometry through these, compared with the 
Johnson- Cousins, filters. However, there have been some 
suggestions of offsets between SPS models and SDSS ob- 
servations. In particular, studies of luminous red galaxies 
(LRGs, which are expected to have a predominantly old, 
metal rich, stellar population) have shown significant off- 
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sets from the model predictions, part i cularly around th e 
^n-bands (e.g. lEisenstein et al.|[200lHWake et al.ll2006D . 
Several expla nations have bee n proposed to rectify this 
offset. Mara ston et al.l ([20091 ) proposed that slight dis- 
crepancies between theoretical and empirical stellar li- 
braries may explain the observed offset. They found 
that models which include a minority population of metal 
poor stars (3% of the total) and that are based on em- 
pirical, rather than theoretical, libraries better represent 
the observed colors a t all redshifts considered. How- 
ever, IConrov fc Cil^ (p310) were unable to explain the 
LRG colors using these different stellar libraries. Instead, 
it has also been suggested that a small population of 
younger or relatively metal poor stars in these galaxies 
could explain the SDSS observations. The addition of 
other stellar populations such as blue/extreme horizon- 
tal branch (HB) stars and blue straggler stars (typically 
excluded in such models) could also help to explain the 
offsets. Some discrepancies have also been suggested be- 
tween SSP models and the (resolved) SDSS color magni- 
tude diagrams o f 17 of the Milky Way's globular clusters 
(|An et all 12001) . By comparing t heir fiducial sequences 
with the theore tical isochrones of iGirardi et al. (200 J), 
lAn et all (|2008l ) demonstrated that the models could not 
be simultaneously fit to both the main sequence and red 
giant branch accurately. 

In this paper, we compare SSP models with the inte- 
grated colors of globular clusters in M31. This cluster 
system is the only large sample of clusters for which ac- 
curate integrated photometry is available from the SDSS 
survey. It therefore provides a unique location in which 
to compare SDSS photometry of star clusters with SSP 
models. In section 2 we describe the (previously pub- 
lished) data available for MSl's clusters. The SPS models 
considered in this paper are introduced in section 3. Sec- 
tion 4 compares these models with the cluster data and 
section 5 considers possible explanations for observed off- 
sets. 

2. M31 GLOBULAR CLUSTER DATA 

M31 hosts the largest cluster population in our local 
group. The galaxy's proximity to us (~780 kpc) means 
its clusters' stellar populations are unresolved in typi- 
cal ground based observations. However, the integrated 
emission of these clusters can be accurately measured 
with relatively short exposures (compared with cluster 
systems around more distant galaxies). In this paper, 
we conside r all globular c l usters listed as 'class 1' in the 
catalog of iPeacock et al.l ()2010r P1. These clusters have 
non-stellar surface brightness profiles and are confirmed 
to be at the distance of M31 based on spectroscopy or 
high resolution HST observations. This dataset is there- 
fore unlikely to contain significant contamination from 
non-cluster sources. This sample also excludes the re- 
cently identified population of young clusters in M31. 
Young clusters are classified as those clusters having 
colors bluer than the old Milky Way glob ular clusters 
(|Fusi Pecci et al.l 120051 : IPeacock et al.''2010^ or based on 
spectroscopy of the clusters (Caldwell ct al. 2009). The 
catalog includes 416 class 1 clusters, located across the 
galaxy 

The SDSS covered a large region of sky in the direc- 
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tion of M31 as part of a supplemental run to the main 
survey. These observations have previously been used to 
provide u and griz-hand photometry for 73% and 92% of 
M31's known globular clusters, respectively (fewer clus- 
ters are detected in the it-band due to the increased ex- 
tinction and lower sensitivity of the survey at these wave- 
lengths). The photometry was performed using Sex- 
TRACTOR. The colors were obtained through 4" aper- 
tures and total magnitudes measured using an aperture 
in the range 2.8-10.6", depending on the size of the clus- 
ter. The photometry was calibrated using the standard 
SDSS pipeline calibration. As such, the colors are on the 
AB magnitude system. The SDSS calibration is thought 
to be in good agreement with this system for the griz- 
bands, although a slight correction has b een proposed 
for th e M-band (such that MAR=Mqn^q-0.04: lBohlin et al.l 
I20pi|). For full d e tails of these data we refer the reader 
to IPeacock et al.l ()2010[ ). 

Spectroscopic metallicities for 200 of thes e clusters 
are a vailable from the collated catalog of iFan et all 
()2008l ) . These meas u remen ts are taken fro ir i the stud- 
ies of iHuchra et al.] (|1991f) . IBarmbv et all (|2000[ ) and 
IPerrett ct al. (2002). Where clusters are present in more 
than one of these studies, prefere nce was given to th e 
larger and more recent catalog of IPerrett et al.l ()2002D . 
iFan et all ()2008l ) also provide estimates for the reddening 
of these clusters. The extinction in each band was cal- 
culated for each cluster using the IFan et al. (2008) red- 
dening values and the extinction curves of Cardclli ct alj 
(,1989i) , as presented in table 6 of iSchlegel et al. (19981)7 
The reddening of individual clusters in M31 is variable 
across the galaxy. Therefore, we only consider those clus- 
ters for which reddening estimates are available. This re- 
sults in a sample to 200 of M31's globular clusters which 
have photometry, spectroscopic metallicities and redden- 
ing estimates. To ensure the cluster colors used are rela- 
tively robust to variations between the true and assumed 
extinction curves along their line of sight, we only con- 
sider relatively low extinction clusters, with E(B-V)<0.2. 
Our final sample is comprised of 140 clusters. 

3. SIMPLE STELLAR POPULATION MODELS 

In this paper, we consider three of the commonly used 
SPS models, as well as some recent updates to these mod- 
els. We summarize these below, but for details of the 
models, we refer the reader to the ci ted papers. 

Firs tly, we consider the models of iBruzual fc CharlotI 
(pool hereafter BCO30. BC03 provide models based 
on three stellar evolutionary prescriptions. Of these op- 
tions we take the models based on the 'Padova 1994' 
tracks. This model combines these tracks with either 
the STELIB or BaSeL 3.1 stehar libraries. The colors 
of these models (as a function of age and metallicity) in 
the SDSS filters are taken from the tables provided by 
BC03. Models based on two different initial mass func- 
tions (IMFs) are available. Only a small difference is ob- 
served between the different IMFs, however, the model 
with a Chabrier IMF provides a slightly better fit to 
our data than the Salpeter IMF. We therefore select the 
model with this IMF. 

We also consider the models of iMarastoil ()2005l 

taken from; |http; / / www2.iap.fr /users /chariot /bc2003 / galaxev] 
_download.html 
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hereafter M05j3- These models are based on the 
theoretical BaSeL stellar atmospheres ()Kuruczl 119791 : 
iLeieune et al.lll998D comb ined with models for M giants 
from iBessell et aLl ()1989t ) and empirical colors for the 
thermally-pulsating asymptotic giant branch (TP-AGB) 
stars. M05 also contains some treatment of the HB star 
component, providing models for stellar populations with 
either red or blue HB morphologies. The ugriz photome- 
try for these models were taken from the tables provided. 
As demonstrated later, the model which contains a blue- 
HB morphology and Kroupa IMF provides the best fit to 
our data and is the primary model used in our analysis. 

In addition to M05, we also consider the updated mod- 
els of Maraston & Stromback (MNRAS submitted,2011). 
These models are based on different stellar libraries to 
those presented in M05. Unfortunately, these libraries 
cover a more limited parameter space than the Kurucz 
library. The m o del b ased on the STELIB library of 
iLe Borgne et al.l (|2003D only includes [Fe/H]>0.01, and 
is only applicable t o me tal rich glo bular cl u sters. The 
MILES libraries of Islnchcz-Blazqu ezlit al.l (I2006D and 
the ELODIE library of Prugniel et al I (|2007D "~provide 
better metallicity coverage but only cover a limited wave- 
length range, meaning that they are primarily applica- 
ble to the SDSS g and r-bands. We consider only the 
updated models based on the MILES and ELODIE li- 
braries (hereafter referred to as MSll-MILES and MSll- 
ELODIE). The photometry for these models, as a func- 
tion of age and metallicity, was tabulated by convolving 
the SDSS filter responses with the model SEDs. 

Finally, we consider models based on the latest Padova 
tr acks and obt ained fro m the w eb interface CMD 2.3, 
at http://stev.oapd. inaf .it/cmd[ We extracted two sets 
of models iiTthe "SDSS ugriz filters. The first model 
is that of iMarigo et al.l (l2008f) based on the tracks of 
iGirardi et al.l (|2000l . I200l . T his includes the bo lometric 
corrections for carbon stars of lLoidl et al.l ()200lD and was 
calculated with no circumstellar dust. We refer to this 
as the Padova model. The second model is based on the 
same tracks, but with a correction applied to the TP- 
AGB tracks b ased on the new empirical calibration of 
IGirardi et al.l |2010,. their case B correction). We refer 
to this as the corrected Padova model. 

4. SSP AND GLOBULAR CLUSTER COLORS 

Figure [1] shows the ugriz colors predicted by the SSP 
models, discussed in section 3, and the intrinsic ugriz 
colors of M31's clusters, as a function of metallicity. 
The gray pluses show the dereddened colors of M31's 
clusters with E(B-V)<0.2. Higher extinction clusters 
show greater scatter in their color-metallicity relation- 
ships. However, no significant offset is observed between 
the intrinsic colors of clusters with high and low extinc- 
tion. The model tracks plotted in this figure are for 12 
Gyr SSPs from BC03 (blue line), M05 with a blue-HB 
morphology (red line), Padova (green dashed line) and 
corrected Padova (green solid line). A model based on 
combining the M05 and MSll-MILES model is also in- 
cluded in this figure (dashed red line). However, discus- 
sion of this model is saved for section 5.3. The BC03 
track is based on a Chabrier IMF while the other models 

taken from: http://www.icg. port. ac.uk/~maraston/Claudia?s] 
.Stellar _Population_Model.Mml 



are based on a Kroupa IMF. 

The most striking feature, in figure [1] is that all of the 
SSP models predict significantly redder colors than ob- 
served in g-r (with A((7-r)^0.1). While all models are 
offset in this color, the corrected Padova model best rep- 
resents the g-r-metallicity relationship, with the other 
models diverging from the data at low metallicity. How- 
ever, this model is also offset from the data. The offset 
in g-r is at similar wavelengths to that previously ob- 
served between the BC03 and M05 models and the B-V 
colors of globular clusters (e.g. Maraston 2005). How- 
ever, the effect appears to be more pronounced in the 
g-r filters. This shift in g-r is also similar to that ob- 
served between these model s and SDSS obs e rvatio ns of 
LRGs (see e .g. figure 1 of [ Maraston et al.l ()2009D and 
figure 11 of iConrov fc Gunn (2010)). Potentially this 
offset could also be produced by an offset in the metal- 
licity measurements of the clusters. However, we note 
that systematically reducing the metallicity of the clus- 
ters to explain this offset tends to produce new offsets in 
the other color-metallicity correlations. Possible expla- 
nations for the offset in g-r are discussed further in the 
next section. 

Better agreement is found between the other colors ob- 
served and those predicted by the SSP models. It can be 
seen that the corrected Padova model represents all of 
the observed colors better than the original model and is 
in reasonable agreement with all colors except g-r and (to 
a lesser extent) r-i. The BC03 and M05 models predict 
redder g-z colors than observed. These models also pre- 
dict trends between r-i and metallicity that are slightly 
shallower than that observed. However, it should be 
noted that the r-i color covers a relatively small range, so 
that the effect, which is small in absolute color offsets, is 
nonetheless obvious in the plot. The u-g colors are gen- 
erally consistent with the model predictions. We note 
that the observed relationship between u-g and metal- 
licity appears to flatten at low metallicities. However, 
the other colors are consistent with a linear relationship 
between the color and metallicity. In particular the u-z 
color shows a strong correlation with metallicity and is 
consistent with the colors predicted from the BC03, M05 
and corrected Padova models. 

The red dashed lines in figure [T] represent the combina- 
tion of the MSll-MILES model (in the g and r-bands) 
and the M05 model (in the u,i and z-bands). As dis- 
cussed in section 5.3, this update to the M05 model is in 
better agreement with the observed g-r and g-z colors. 

5. MODEL-DATA DISCREPANCIES 

Comparisons between the observed and predicted col- 
ors of clusters provide a useful diagnostic with which 
to compare different models and highlight potential im- 
provements. Here we consider some of the possible expla- 
nations/remedies for the SSP models predicting redder 
colors than observed around the g and r-bands. 

5.1. Reliability of cluster colors 

Firstl y, we consider the reliability of M31's cluster 
colors. IPeacock et al.l ()2010f) have previously compared 
these data with independent photometry of M31's clus- 
ters obtained in the Johnson-Cousins filters. This was 
do ne using the photo metric transformations presented 
in iJester et al.l (|2005D . Using these transformations, it 
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Fig. 1. — Metallicity-color relationships of M31's globular clusters. Only low extinction clusters, with E(B-V)<0.2, are plotted (gray 
pluses). Assuming a[CardcUi ct al. (1989) extinction curve, this corresponds to E(m-3)<0.27, E(sr-r)<0.21, E{r-i)<0.13, E(j-z)<0.12, E(g- 
z)<0.46, E(m-z:)<0.73. Lines show the predicted colors as a function of metallicity for 12 Gyr SSP models from BC03 (blue), M05 (red), 
MS11-MILES+M05 (red dashed), Padova (green dashed) and Padova-corrected (green). 



was shown that the ugriz colors a re consistent with th e 
previous UBVRcIc photometry of IBarmbv et all (|2000D . 
A shght offset was found in the g-r color (in the sense 
that (?-r=g-rBarmby+0.03). However, this offset is within 
the quoted uncertainty of the transformations used to 
predict the g-r color from the Barmby data. Therefore 
the two datasets are consistent, to within the accuracy 
of the transformations. We also note that the slight shift 
in g-r is smaller than, and in the opposite direction to, 
that observed between M31's cluster colors and the colors 
predicted from SSP models. 

It is also desirable to compare the colors of M31's clus- 
ters with the colors of the Galactic globular clusters. Un- 
fortunately, integrated photometry of the Galactic clus- 
ters is not currently available in the ugriz bands. Instead, 
the ugriz colors of these clusters can be estimated by 
transforming photometry in the UBVRcIc bands (taken 
from the 'December 2010' version of the Harris catalog; 
[Harris 1996 ) to the ugri z ban ds using the transforma- 
tions of iRodgers et al.l (j2006( ) . These transformations 



are based on comparison between standard stars that are 
common to the UBVRcIc and u'g'r'i'z' filter systems. A 
comparison between the colors of the Milky Way's and 
M31's cluster systems is shown in figure [5] In all col- 
ors, fewer red (metal rich) clusters are observed in the 
Milky Way sample. This is likely due to only including 
Milky Way clusters with E(B-V)<0.4. While this cut 
ensures the cluster colors are relatively reliable, it also 
excludes the majority of the metal rich (redder) clusters 
in the Galaxy. The g-r and r-i colors are found to be 
similar between the two cluster systems. It can be seen 
that M31's clusters are found to extend to bluer u-g col- 
ors than the Milky Way's. However, the transformation 
between the ug and UB filters is relatively complicated 
because these filters cover the Balmer break region of 
the stellar spectra. This produces a jump in the trans- 
fo rmation between the u-g and U-B colors (see figure 3 
of IRodgers et al.l I2006D and is the likely reason for the 
transformed Milky Way cluster u-g colors grouping in 
a relatively narrow range. It should be noted that the 
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Fig. 2. — Colors of clusters in M31 (black points), transformed colors of Milky Way clusters with E(B-V)<0.4 (gray/black crosses) and 
g-z color of clusters in the VCS (gray points). Absolute g-band magnitudes assume a distance modulus of 24.4 and 31.1 for M31 and the 
Virgo cluster, respectively. 



transformations of IRodgers et al.l (|2OO60 are based pri- 
marily on main sequence stars. The integrated emission 
from these globular clusters will be slightly different to 
this, including a significant component from cooler red 
giant stars. It is therefore likely that the difference in u- 
g is due to the photometric transformations used, rather 
than suggesting an intrinsic difference between the colors 
of the two cluster systems. 

Relatively little data for other cluster systems has been 
published in the ugriz filters. However, the g-z color of a 
large sample of clusters is available from t he HST/ACS 
Virgo Cluster Survey (VCS; as published bv lJordan et al.l 
120091 and obtained from the VizieR archive). The VCS 
obtained photometry through the F475W (SDSS g) and 
F850LP (SDSS z) fihers of globular clusters in 100 dif- 
ferent Virgo cluster galaxies. In the bottom right panel 
of figure [2] we compare the colors of clusters in the VCS 
(smaU gray points) with M31's clusters (black points). It 
can be seen that M31's clusters span a similar range of 
g-z colors. More red (metal rich) clusters are observed 
in the VCS dataset, compared with M31. However, this 
is unsurprising given the large number of clusters in el- 
liptical galaxies in the VCS dataset (these galaxies are 
known to host relatively large populations of metal rich 
clusters) . 

A potential source of uncertainty on the intrinsic colors 
of MSl's clusters is reddening, which varies significantly 
due to M31 itself. These data were dereddened using the 



values determined bv lFan et al.l (|2008[ ). However, during 
the course of this work, a new spectros copic survey o f 
M31's clusters was published by Caldwei reTall ()2011l) . 
This provides independent measurements for the metal- 
licities and redden ings of these clust ers. Comparing these 
data with that of iFan et all ()2008f ). we note only small 
offsets between the two data sets of A[Fe/H]~0.1 and 
AE(B-V)^0.05. A dopting the met a llicity and reddening 
values presented bv lCaldwell et al.l ()2011l ) actually tends 
to increase the observed offset between the model and 
data colors. It is therefore unlikely that the offset (seen 
in figure [1]) is due to systematic errors in either the red- 
dening correction or spectroscopic metallicities. In this 
work, we have kept the reddening values presented in 
I Fan et all (2008) because they are found to produce less 
scatter in both color-metallicity and color-color plots. 

Considering these comparisons, it appears unlikely 
that the shifts we observe between the data and mod- 
els are due to the calibration of the cluster colors. 



5.2. Additional stellar populations 

A potential explanation for offsets between the ob- 
served and predicted colors of these clusters is a dif- 
ference in the underlying stellar populations. Certain 
phases of stellar evolution are not included in these SPS 
models, or are known to be relatively poorly treated. One 
such population are the HB stars, which are notoriously 
difficult to model. In figure [3l we show the M05 models 
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with red-HB star (dashed red hnes) and blue-HB star 
(sohd red hnes) morphologies. It can be seen that the 
addition of blue-HB stars does produce bluer integrated 
colors and therefore represent the data better, partic- 
ularly at low metallicity (where clusters are known to 
have bluer- HB star populations. iSandage fc Wallersteinl 
IT960t) . However, even for the metal poor clusters, the 
shift is far less than that required to explain the observed 
colors. Including more blue/extreme-HB stars could po- 
tentially push the models to even bluer colors, but it is 
unlikely this could produce the (large) observed offset. 
Blue straggler (BS) stars are also known to be present 
in globular clusters and are generally excluded from SSP 
models. These bright, blue, stars would tend to shift the 
models towards bluer colors. However, clusters typically 
contain only a small fraction of blue stragglers, so it is 
unlikely they will produce the large offset seen in most of 
M31's clusters. For example, if we assume that a cluster 
contains approxima tely 1 blue straggler star per 2000Mo 
(jKnigge et al.ll2009l) and that a typical blue straggler has 
g^2A and g-r~0 |Leigh et al.l [20071. we would expect 
the addition of such a population to make the model 
of an old, metal poor cluster bluer by only g-r^O.Ol. 
Therefore, while the addition of these stellar populations 
can produce bluer cluster colors, it appears unlikely that 
they can produce the large offset observed in M31's clus- 
ters. The addition of these hot populations would also 
be expected to influence the u-g color of these clusters. 
However, relatively good agreement is found between the 
observed and expected u-g colors, compared with the g-r 
colors (see figure [T|). 

It has also been proposed that some of the observed 
offsets between SPS models and LRG observations, can 
be explained by adding minority populations of metal 
poor stars ( Maras tpn et al. 209% Conr oy & Gunn 2010) 
or young stars (jConrov fc'Gunni (201 Ol) to the expected 
old, metal rich, stellar population. While such popula- 
tions may be possible in these galaxies, age and metal- 
licity spreads on this scale are not observed in globular 
clusters. It is therefore very unlikely that the observed 
colors of M31's clusters are due to minority populations 
of younger or lower metallicity stars. 

5.3. Stellar libraries 

Another important factor in comparing SPS models 
with observations is the stellar libraries that are used to 
convert the models in to observable spectra. Previously, 
IMaraston et al.l (j2009D have proposed that the observed 
offset between the M05 models and observations of LRGs 
in the SDSS, could be explained by using updated mod- 
els based on empirical, rather than the theoretical, stellar 
libraries. By using these empirical libraries, and includ- 
ing a minority population of low metallicity stars, they 
found good agreement with the ^ri-band observations of 
these galaxies at all redshifts cons i dered . However, in a 
subsequent paper [Conrov fc GunnI (|2010[ ) were unable to 
explain this offset by using the same empirical libraries 
in their models. By using these M31 data, we can com- 
pare models based on different stellar libraries with the 
relatively simple stellar populations of globular clusters, 
over a greater spread in metallicity. 

In figure[3]we compare the g-r colors of M31's clusters 
with the original SSP models of M05, based on the theo- 
retical Kurucz stellar library (red lines) . Considering the 




[Fe/H] 

Fig. 3. — Comparison between models using different stellar li- 
braries. As in figure [T] M31's clusters are shown in gray. The 
black line represents the average color of MSl's clusters (based on 
grouping clusters in to 15 metallicity bins). The red lines show 
the original models of M05 for 12 Gyr SSPs with blue and red- 
HB morphologies (solid and dashed lines, respectively). The upper 
dashed and upper solid red lines are models based on a Salpeter 
IMF, while the lower red lines are based on a Kroupa IMF. The 
green line shows the MSll-ELODIE model for an 11 Gyr SSP. The 
solid blue line shows the MSll-MILES model f or a 12 Gyr SSP and 
the dashed blue line the lVazdekis et al.l | |2010I '| model. 

different versions of the M05 models, it can be seen that 
the models which include a blue-HB morphology (solid 
red lines) provide a better fit to the observations. The 
M05 models based on a Kroupa IMF (lower red dashed/ 
solid lines) produce bluer colors than those based on a 
Salpeter IMF and are therefore in better agreement with 
our observations. However, all of the M05 models pro- 
duce redder colors than observed, for all metallicities. 

The solid blue line in figure |3] shows the colors pre- 
dicted from the MSll-MILES model. This model is 
based on the same evolutionary code as the M05 models 
but uses the empirical MILES stellar library. It can be 
seen that this model is in relatively good agreement with 
the observed g-r colors for all metallicities. This is con- 
sistent with the findings of IMaraston et ahl (j2009[ ). who 
found that these models also provide a good fit to the 
SDSS LRG observations. Compared with the theoretical 
libraries, the MILES stellar library has relatively lim- 
ited wavelength coverage. Therefore, the MSll-MILES 
model can not be directly applied to the other SDSS 
bands studied. However, we note that if these shifts in 
the g and r-bands are applied to the original M05 mod- 
els, it also improves the fit to the observed g-z colors 
of M31's clusters. This updated model is represented in 
figure [T] as the dashed red line. Reasonable fits are also 
obtained to the u-g color, although the updated models 
appear to be slightly too red compared with the metal 
poor clusters. The trend between r-i and [Fe/H] remains 
slightly flat compared with the data. Again we note that 
r-i only covers a relatively small range and that the off- 
set between the data and models is small. However, it is 
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possible that the implied correction to the model spec- 
trum around the g and r-bands may also have a (smaller) 
influence on the i and w-bands. 

An independent stellar population model, that is also 
ba sed on the MILES ste llar library, has been presented 
by iVazdekis et aP (j2010l ). The g-r color of this model 
was obtained using the data access tooQ. The blue- 
dashed line in figure [3] shows the colors predicted from 
this model. It can be seen that, compared with the M05 
models, this model provides a better fit to the metal 
poor clusters. However, it is still redder than the ob- 
served c olors, especially at higher metallicities. Given 
that the IVazdekis et al.i (|2010D and MSll-MILES mod- 
els are based on the same stellar library, it is interest- 
ing that they predict significantly different colors around 
solar metallicity. The offset between the two models 
demonstrates the complexity of predicting integrated col- 
ors from SPS models, which are highly dependent upon 
the evolutionary prescriptions used and the implemen- 
tation of the stellar libraries, as well as the libraries 
themselves. Our comparison suggests that the use of the 
MILES library can help to explain the offset between 
the Maraston models and observations. However, us- 
ing this library does not appear to explain the similar 
offset observed betwe en o bservations and t h e mo dels of 
IVazdekis et all ([201(11 ) and lConrov fc GunnI (poTfih . The 
reason for the observed offset between these models and 
observations remains unclear. The different predictions 
of the MSll-MILES and Vazdekis/Conroy models also 
remains an open issue for the modeling community. How- 
ever, it may be related to differences in the treatments of 
the red giant branch/ asymptotic giant branch stars in 
the Maraston and the (Padova based) Vazdekis/Conroy 
models. 

Figure 121 also shows the MSll-ELODIE model. This 
model is similar to the M05 and MSll-MILES models, 
but is based on the ELODIE stellar library (green line). 
It can be seen that this model also provides a better 
fit to the low metallicity clusters. However, it still pre- 
dicts redder colors than those observed for the higher 
metallicity clusters. The difference between this model 
and the MSll-MILES model must be due to differences 
between the stellar libraries. A likely reason for this is 
difficulties in fiux calibrating the stellar spectra. Indeed, 
MS 11 demonstrate that significant differences are found 
between the spectra of stars that are present in both 
the MILES and ELODIE libraries. The variation in the 
colors predicted by the otherwise similar M05, MSll- 
MILES and MSll-ELODIE models highlights the impor- 
tance of stellar libraries to the predictions of SPS models. 

To highlight the differences in these models, we show in 
figure m the spectra predicted for stellar populations with 
solar metallicity. For clarity, we have binned the higher 
spectral resolution models to that of the M05 models 
and scaled all fluxes to that at 5500A (close to the edge 
of the g and r-band filters). This figure demonstrates 
that the MSll-MILES model produces slightly higher 
fiuxes than the M05 model at most wavelengths covered 
by the SDSS g-h&nd filter. Conversely, this model pro- 
duces slightly lower fluxes across the SDSS r-band fllter. 
This leads to the bluer g-r color predicted by this model. 

* provided at http://milcs.iac.es/ 
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Fig. 4.— Spectra of 12 Gyr SSPs with [Fe/H]=-f 0.00 produced 
from the models of M05 (red), Vasdekis (black), MSll-MILES 
(blue) and MSll-ELODIE (green). All models have been scaled 
by their flux at 5500A, this is close to the intersection between the 
g and r-band filters. The models have also been binned to the 
spatial resolution of the M05 models. For comparison, the SDSS 
filter transmission curves are shown in the top panel. 

In comparison to the MSll-MILES model, the use of 
the ELODIE library in the MSll-ELODIE model pro- 
duces slightly less flux in the i^-band filter. However, the 
biggest difference between these models is around 6500A, 
where the MSll-ELODIE model produces higher fiuxes. 
In the bottom panel of figure [H it can be seen that, 
the MSll-MILES and Vazdckis models produce similar 
fiuxes across the g-hai\A. However, the Vazdekis model 
predicts increasingly higher fiuxes at longer wavelengths. 
This produces the the redder g-r color predicted by the 
Vazdekis model, compared with the MSll-MILES model. 

6. CONCLUSIONS 

We have demonstrated a significant offset between the 
colors of M31's globular clusters and those predicted by 
some SSP models. The biggest discrepancy is found 
around the SDSS g and r-bands. The offset observed 
in g-r is similar to that previously observed in Sloan 
LRG datasets, but is confirmed here to also effect the 
relatively simple stellar populations of globular clusters, 
across a range of metallicities. Identifying this offset in 
globular clusters suggests that it is unlikely to be due to 
a small fraction of young stars, as has previously been 
proposed in the LRGs. Including other minority popu- 
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lations such as cxtrcmc/bluc-HB stars or blue straggler 
stars can produce clusters with bluer colors. However, 
such populations would not be expected to produce the 
large observed offset in only g-r. 

We have shown that the g-r colors predicted by the 
MSll-MILES model are consistent the observed clus- 
ter colors, for all metallicities. However, the same 
models based on the theoretical, Kurucz, and empiri- 
cal, ELODIE, libraries predict significantly redder colors 
than observed. This suggests that the observed offset 
in g-r can be explained by differences between the stel- 
lar libraries used. It should be noted that the Vazdekis 
model, also based on the MILES library, is still offset 
in g-r. The reasons for the difference between the MSll- 
MILES and Vazdekis models are currently uncertain and 
worthy of future investigation. However, we note that the 
MSll-MILES model is the only model-library combina- 



tion considered that gives a reasonable representation of 
the observed colors of M31's clusters. 

Considering such shifts between models and observa- 
tions provides an important tool for refining SPS mod- 
els. Such offsets should be considered when inferring the 
properties of globular clusters and galaxies from their 
integrated emission. 
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